Thermally activated energy and flux flow Hall eff'ect of Fei+y(Tei_j:Sx)z 
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Thermally activated flux flow (TAFF) and flux flow Hall efl'ect (FFHE) of Fe(Te,S) single crystal 
in the mixed state are studied in magnetic fields up to 35 T. Thermally activated energy (TAE) 
is analyzed using conventional Arrhenius relation and modified TAFF theory which is closer to 
experimental results. The results indicate that there is a crossover from single-vortex pinning region 
to collective creep pinning region with increasing magnetic field. The temperature dependence of 
TAE is different for H / jah and H/ /c. On the other hand, the analysis of FFHE in the mixed state 
indicates that there is no Hall sign reversal. We also observe scaling behavior \pxy{H)\ — Apxx{H)^ . 

PACS numbers: 74.25.Wx, 74.25.F-, 74.25. Op, 74.70.Dd 
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I. INTRODUCTION 

The discovery of iron based materials has generated 
enormous interests in the field of superconductivityii~— 
Due to similar layered structure to cuprate oxides and 
rather high Tc, iron based superconductors could host 
rich vortex phenomena in the mixed statej^"— Recent 
vifork suggests that vortex properties of iron based mate- 
rials seem to be similar to the cuprate superconductors 
since magnetic flux coUective pinning and creep region 
as well as fishtail effects (second peak effect) have been 
observed in LnOFeAs (Ln=rare earth elements, 1111 
system) and AFe2As2 (A~alkaline earth elements, 122 
system) 

Iron based superconductors, FeSei-a;, Fei+yTei-xSe^:, 
and Fei+yTei_:rSa; (11 system))^J^~— are of interest 
both for the technological applications and for the un- 
derstanding of the vortex properties in the mixed state 
due to rather simple structure and nearly isotropic upper 
critical fieldi^~— Only a limited amount of information 
on the vortex behavior in single crystals of 11 system is 
available until now, mainly focusing on thermally acti- 
vated flux flow (TAFF) region in Fe(Te,Se),25^ Normal 
carriers in the vortex core, which experience a Lorentz 
force, can lead to normal Hall effect in mixed state. On 
the other hand, flux flow can also induce Hall effect in 
the mixed state. In detail, when applying a transport 
current, the flux lines will experience the Lorentz force, 
F = ij X B, where j is the supercurrent density and B 
is the magnitude of magnetic induction. The motion of 
magnetic flux lines produces a macroscopic electric field 
E which is given by E = — -^v x B, where v is the velocity 
of vortex motionj^ The vortex motion along the Lorentz 
force (perpendicular to j) gives the dissipative field (E||j) 
and generates the flux flow resistivity, whereas the vor- 
tex motion along the direction of supercurrent results in 
the HaU electric field (E±j). Therefore the flux flow HaU 
effect (FFHE) is a sensitive method to study the vortex 
dynamics. 

There are two exotic phenomena in connection with 
FFHE in cuprate superconductors. One is a sign reversal 



of the Hall resistivity pxy{H) below Tc.= This anomaly 
has also been observed in some conventional supercon- 
ductors, e.g. amorphous MoSis,^^ and 2H-NbSc2.^2 The 
sign change is not expected within the classical Bardeen- 
Stephen-^ and Nozieres-Vinen'^^ theories of vortex mo- 
tion, which predict that the Hall sign in the supercon- 
ducting and normal state should be the same. Sev- 
eral models have been proposed for interpreting this 
anomaly;^"— however its origin remains a controversy. 
Another phenomenon is a scaling law between pxy{H) 
and the longitudinal resistivity Pxx[H) in the supercon- 
ducting transition region, i.e., \pxy{H)\ = Apxx{H)^ with 
different values of /3 for different materials i^^^ 

In this paper, we study the vortex properties of 
Fei.i4(i)(Tco.9i(2)So.o9(2))^ single crystal via TAFF resis- 
tivity and FFHE in the mixed state. The temperature 
dependence of TAE is different for H\\ab and H\\c. Fur- 
thermore, there is a crossover from single- vortex pinning 
region to collective creep pinning region with increasing 
magnetic field. On the other hand, there is no sign reveral 
and we observe scaling behavior for FFHE. 



II. EXPERIMENT 

Single crystals of Fe(Te,S) were grown by self flux 
method and their crystal structure was analyzed in 
the previous report.^- The elemental analysis of the 
crystal used in this study showed the stoichiometry is 
Fe:Te:S=1.14(l):0.91(l):0.09(2) and we denote it as S-09 
in the following for brevity. Electrical transport measure- 
ments were performed using a four-probe configuration 
with current flowing in the ab-plane of tetragonal struc- 
ture in dc magnetic flelds up to 9 T in a Quantum Design 
PPMS-9 from 1.9 to 200 K and up to 35 T in a He3 cryo- 
stat system with resistive magnet down to 0.3 K at the 
National High Magnetic Field Laboratory (NHMFL) in 
Tallahassee, FL. Hall contacts with typical misalignment 
of less than 0.1 mm were used. At each point the Hall 
voltage was measured for two directions of the magnetic 
field, which is always perpendicular to current direction. 
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III. RESULTS AND DISCUSSIONS 

Fig. l(a,b) show the resistivity p{T,H) of S-09 near 
the superconducting transition region for H \\ ab plane 
and if II c axis. With increasing magnetic fields, the re- 
sistivity transition widths are broadened gradually. The 
onset of superconductivity shifts to lower temperatures 
for both magnetic field directions, but the trend is more 
obvious for H\\c than H\\ab. According to the TAFF 
theory, the resistivity in TAFF region can be expressed 

p = {2:yQLB/J)exp{-JcoBVL/T)sinh{JBVL/T) (1) 

where vq is an attempt frequency for a flux bundle 
hopping, L is the hopping distance, B is the magnetic 
induction, J is the applied current density, Jco is the 
critical current density in the absence of flux creep, V 
is the bundle volume and T is the temperature. If J is 
small enough and JBVL/T <C 1, we obtain 

p = {2pJJ /T)exp{-U /T) = pofexp{-U/T) (2) 

where U = JcqBVL is the thermally activated energy 
(TAE) and pc = vqLB / Jco, which is usually considered 
to be temperature-independent. For cuprate supercon- 
ductors, the prefactor 2pcUlT is usually assumed as a 
constant po/r therefore, lnp{T, H) = lnpQf — U{T, H)/T, 
where H is the external magnetic field. On the other 
hand, according to the condensation model, ■^^ U (T, H) = 
where He is the thermal critical field, ^ is 
the coherence length, t = T/Tc [Tc is the superconduct- 
ing transition temperature), and n depends on the di- 
mensionality of the vortex system with the range from 
to 3. Since iJc oc 1 — and ^ oc (1 — t)~^/'^ near 
it is obtained that U{T, H) = Uo{H){l -tY , where 
g = 2 — n/2. It is generally assumed that U{T,H) = 
Uo{H){l — t), i.e. n = 2, and the Inp — 1/T becomes Ar- 
rhcnius relation, lnp{T,H) — lnpo{H) — Uo{H)/T, where 
lnpo{H) — Inpof + Uo{H)/Tc and Uo{H) is the apparent 
activation energy. Furthermore, it can be concluded that 
-dlnp{T,H)/dT-^ = Uo{H). Hence, the Inp vs. 1/T 
should be linear in TAFF region. The slope is Uo{H) 
and its y- intercept represents lnpo{H). 

In the Figs. l(a,b), the solid lines show the Arrhenius 
relation in TAFF region. Note that the results are shown 
in the common logarithmic scale in the figures, but we 
calculate them in the natural one. All linear fits intersect 
at approximately the same point Tcross, which is about 
8.63 and 8.34 K for H \\ ab and H \\ c, respectively. 
Assuming the temperature dependences of pab(T,H) at 
two different magnetic fields {Hi and H2) can be fitted 
by Arrhenius relation, according to above discussion, we 
get lnp{T,Hi) = Inpof + Uq{Hi)/T, - Uo{Hi)/T and 
lnp{T,H2) = Inpof + Uo{H2)/T, ~ Uo{H2)/T. When 
p{T,Hi) = p{T,H2), it can be obtained T = there- 
fore, ideally, all the lines at different fields should be 
crossed into same point, Tcross^ which is equal to Tc- 
According to the conventional analysis, the Inpof and Tc 
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FIG. 1. (a) and (b) Longitudinal resistivities p{T,H) of S-09 
in different magnetic field directions for _ff||a6 and H\\c below 
lOK, respectively. The corresponding solid and black dashed 
lines are fitting results from the Arrhenius relation and eq. 
(3). (c) Inpo vs Uo derived from Arrhenius relation for _H"||a6 
and H\\c. 



can be obtained from linear fits lnpa{H) and Uq{H) using 
lnpo{H) ^ lnpof + Uo{H)/Tc (shown in Fig. 1(c)). From 
the fitting results, values of pof and Tc are 27.22 ± 4.77 
mQ-cm, 8.72±0.82Kand5.23±I.42 mn-cm, 8.32±0.48 
K for _ff II ab and H \\ c, respectively. The Tc is consis- 
tent with the values of Tcross in the range of errors. It 
seems that the p{T, H) can be fitted with straight lines 
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well. However, close inspection shows that there is rather 
large fitting errors, especially for H \\ ah. The origin of 
the large errors is the Arrhenius relation that can only be 
satisfied in the limited region and this region is narrower 
for II ah. The effects of prefactor and non-linear rela- 
tion of U{T, H) lead to p{T, H) deviating from Arrhenius 
relation (vide infra). 

Fig. 2(a,b) shows the temperature dependence of 
—dlnp{T, H)/dT~^ for both field directions. Because the 
assumptions U{T,H) = Uq{H){1 — t) and pof = const 
lead to -dlnp{T,H)/dT-^ = Uo{H), Uq{H) should be 
a set of horizontal lines. We present this in Fig. 2(a,b) 
over a limited length. Each length covers the temper- 
ature interval used for estimating Uq{H) in the Arrhe- 
nius relation. It can be seen that —dlnp{T, H)/dT~^ 
increases sharply with decreasing temperature, which 
was also observed in Bi-2212 thin films. The cen- 
ter of each Uo{H) horizontal line approximately inter- 
sects —dlnp{T, H)/dT^^ curve and the overlapping re- 
gion is increasing with temperature decrease. This 
shows that each Uo{H) is only the average value of 
its —dlnp{T, H)/dT^^ in the fitting temperature re- 
gion. Hence, the TAE determined from the conventional 
method does not reflect the true evolution of U{T, H) 
with the temperature, particularly for iJ || ah. This con- 
tradiction originates from two basic assumptions intro- 
duced for Arrhenius relation: one is the constant pref- 
actor pof — '^PcU/T and another is the linear relation 
U{T,H) = Uq{H){1 - t). Zhang et al^"^ suggested that 
the temperature-dependent prefactor and nonlinear rela- 
tion of U{T, H)—T should be considered. In the following 
section, we will analyze the resistivity results using this 
more general method. ^S. 

Using the relation U{T,H) = L/o(i?)(l - from eq. 
(2) it can be derived that 




123456789 10 
T(K) 



FIG. 2. Experimental -dlnp/dT'^ data in TAFF region for 
H\\ab (a) and -ff||c (b), respectively. The red solid horizontal 
lines correspond to obtained Uo{H) from Arrhenius relation 
and the blue dashed lines are plotted using eq. (4). The 
parameters are determined via fitting eq. (3) to corresponding 
experimental resistivity data shown in Fig.l (a) and (b). 



Inp = ln{2pJJo) + qln(l - t) - \nT - Uo{l - tf /T (3) 
and 

-dlnp/dT-^^[Ua{\-tf-T][l + qt/{l~t)] (4) 

where pc and Uq are temperature independent and 
derived from Arrhenius relation is used for fitting. There- 
fore, there are three free parameters, q, pc and Uq in 
eq. (3). The fitting results are shown in Fig. 1. From 
Fig. 1 (a,b) , it can be seen that all fits are in good agree- 
ment with experimental data and the results are better 
than Arrhenius relation. This is more pronounced for 
H II ah. Fig. 2(a,b) clearly shows the advantage of eq. 
(3) over Arrhenius relation. The TAFF formula (eq. (1)) 
can effectively capture the upturn trend of —dlnp/dT^^ 
with decreasing temperature when the (linear or nonlin- 
ear) correlations between prefactor as well as U (T, H) 
and T are considered. In detail, when T U (corre- 
sponding to T ^ Uq{\ — tY)^ it can be derived that 
-dlnp/dT-^ = Uo{l-t)i[l + qt/{l-t)] and when g = 1, 
—dlnp/dT^^ = L/qj i-G. Arrhenius relation. Because the 



obtained Uq of S-09 is much smaller than that of cuprates 
superconductors (shown in Fig. 3) and it is comparable 
with temperature, the assumption T <ti U can not be 
satisfied and the temperature dependence of prefactor 
should be considered. 

Fig. 3 presents the Uo{H) and q{H) obtained from 
experimental data fits using eq. (3) at different fields. 
The Uo{H) shows a power law {Uo{H) - iJ"") field de- 
pendence for both directions. For H\\ah, a — 0.12 ± 0.02 
for poH < 5T and a = 1.70 ± 0.30 for hqH > 5T; For 
H\\c, a ^ 0.21 ± 0.03 for hqH < 5T and a = 1.34 ± 0.16 
for pqH > 5T. The weak field dependence of Uq{H) 
in low field for both orientations suggests that single- 
vortex pinning dominates in this region.'' The vortex 
spacing becomes significantly smaller than penetration 
depth in higher fields and we expect a crossover to a 
collective-pinning regime where the activation energy be- 
comes strongly dependent on the field, i.e., the collective 
creep dominance.— Because a is larger than 1 at this 
regime, it is possible that the fiux lines are pinned by the 
collective point defects in the high field region.'*'* Sim- 
ilar crossover has been observed in Nd(0,F)FeAs single 
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FIG. 3. (a) Uo and (b) q as functions of magnetic fields ob- 
tained from fitting the resistivity in TAFF region using eq. 
(3). The opened and filled squares represent Uq for H\\c and 
H\\ab, respectively, while the opened and filled circles show 
corresponding q, respectively. The solid lines in (a) are power- 
law fitting using Uo{H) ~ H~". 



crystalM- The values of q change from about 1 for H\\c to 
2 for H\\ab, independent on the intensity of field for both 
directions. The value oi q = 2 has also been observed in 
many cuprates superconductors.—'^'^ 

It should be noted that although the Fe(Te,S) sample 
has low volume fraction, it should not have obvious ef- 
fect on the analysis of TAFF. Due to the inhomogeneity 
of sample, the conductivity of sample can be expressed 
as a = CTsc + (^normal, whcrc (T sc IS the conductivity in 
superconducting state and (Jnormai is the conductivity 
of normal state. When T < Tc^zero, o'sc is infinite and 
thus cr = (Tsc = oo, i.e., p = and the normal state 
part of sample is short-circuited. On the other hand, the 
resistivity in the TAFF regime is about one to three or- 
ders of magnitude less than normal state2.i^. It means 
that although the conductivity <Tsc in this range is finite, 
it is still much larger than UnormaU and we can obtain 

The Hall effect in the mixed state gives important in- 
sight in the flux flow. In the following section, we will dis- 
cuss the vortex dynamics of S-09 in flux flow region. The 
field dependence of the longitudinal resistivity Pxx{H) 
for H/ /c is shown in Fig. 4(a). Superconductivity is 
suppressed by increasing magnetic field up to 35 T and 
the transition of Pxx{H) are shifted to lower magnetic 
fields at higher temperature. At the lowest measuring 
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FIG. 4. Field dependence of (a) Longitudinal resistivity 
Pxx{H) and (b) Hall resistivity pxy{H) at various tempera- 
tures in dc magnetic fields up to 35T for H\\c. 



temperature (T = O.S-ftT), normal state is recovered from 
superconducting state when field is up to 30 T. Fig. 4(b) 
shows the Hall resistivity at T ^ lOK . It can be seen that 
Pxy{H) = in low field when temperature is below Tc- 
At higher field region close to the superconducting tran- 
sition, the absolute values of Hall resistivity increase and 
gradually reach the Pxy{H) curve obtained in the nor- 
mal state at temperatures slightly higher than T^. The 
Pxy{H) in the mixed state shifts with increasing temper- 
ature to lower fields. All of these features are similar to 
the Pxx{H) results. On the other hand, the normal state 
Pxy{H) curves are very close to linear, except for low 
field parts where there is slight nonlinearity with nega- 
tive curvatures that can be ascribed to skew scattering 
due to excess Fe.*^ 

The sign of the Hall resistivity is positive in the mixed 
state as well as in the normal state indicating hole type 
carriers. There is no sign reversal for pxy{H) in the 
mixed state, which is a typical behavior for hole- and 
electron- type cuprate superconductors below Tj,-^^*'"*^ Be- 
cause Hall conductivity axy{H) [~ Pxy{H)/ {pxx{H)'^ + 

Pxy{Hf) = Pxy{H) / Pxx{H? , whcn Pxx{H) > \pxy{H)\] 

is usually insensitive to disorder by a general argument of 
the vortex dynamics, it is convenient to discuss the Hall 
results using axy{H)i^'^ There are two contributions to 
the axy{H) in the mixed state: 

<Txy{H) = (Txy.niH) + (Jxy.sc{H) (5) 
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where (Jxy,n{H) is the conductivity of normal quasipar- 
ticles that experience a Lorentz force inside and around 
the vortex core. This term has the same sign as the 
normal state and is proportional to H. The second term 
'^xy,sc{H) is an anomalous contribution due to the mo- 
tion of vortices parallel to the electrical current density 
j- 

From the theory based on the time-dependent 
Ginzburg-Landau (TDGL) equation, (Txy,sc{H) cc 1/H 
and it could have a sign opposite to that of (Jxy,n{H)^'-^ 
Furthermore, the <Jxy^sc{H) is the dominant term at low 
field but at higher field <Jxy,n{H) are important and could 
dominate over axy,sc{H). Therefore, if <Jxy.sc{H) has a 
different sign when compared to axy,n{H), it is possi- 
ble to observe a sign reversal in the Hall effect in the 
superconducting state^^i^ as for example in YBCOj^ 
On the other hand, it can be easily seen in Fig. 5(a) 
that the Hall conductivity decreases with increasing field 
and the field dependence of Uxy [H) changes more rapidly 
than 1/H. This suggests that axy{H) is not indepen- 
dent of disorder in the strong pinning regime. Accord- 
ing to the theory proposed by Fukuyama, Ebisawa, and 
Tsuzuki (FET),^^ the sign of cTxy,sc{H) is given by the 
sign of sgn{e)dN{^)/d^\fj,=EF, where sgn[e) is the sign 
of the carrier, A^(/i) is the density of states, is the 
chemical potential, and Ep is the Fermi energy. On the 
other hand, in the phenomenological theory based on 
Ginsburg-Landau equation and its gauge invariancej^ 
the sign of the axy,sc{H) is determined by the signs of 
sgn{e)d\nTc/d^i. In any case, the sign of the Hall ef- 
fect in the mixed state depends on the details of the 
band structure. For a complicated Fermi surface, the 
signs of (Jxy,sc{H) may be different from that of <Jxy,n{H). 
Therefore, contrary to cuprate superconductors, the dif- 
ference of Fermi surface may be one origin of absence of 
sign reversal in Fe(Te,S). On the other hand, cuprate 
superconductors are d-wave superconductors, whereas, 
for Fe(Te,S), the gap function is unknown, but is most 
likely s-wave. This difference could be another origin 
of different contribution of the vortex cores to the Hall 
conductivity!^ 

At high magnetic fields, we observe \pxy{H)\ — 
Apxx{HY scaling (Fig. 4(b)). The values of (3 are in the 
range of 0.9 — 1.0 and increase slightly with temperature. 
A phenomenological model considering the effect of pin- 
ning on the Hall resistivity proposed by Vinokur et al.f^ 
gives the scaling index /3 = 2 which, as well as Hall con- 
ductivity, is independent on the degree of disorder. This 
is believed to be a general feature of any vortex state 
with disorder dominated dynamics. On the other hand, 
based on the normal core model proposed by Bardeen 
and Stephen^ Wang, Dong, and Ting (WDT)^^ de- 
veloped a theory for the Hall effect that includes both 
pinning and thermal fluctuations. In the WDT the- 
ory the scaling behavior is explained by taking into ac- 
count the backflow current of vortices due to pinning. 
Thereby /3 changes from 2 to 1.5 as the pinning strength 
increases)^ This has been observed in irradiated YBCO 
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FIG. 5. (a) Field dependence of absolute values of the Hall 
conductivity \a:cy{H)\ measured at various temperatures in 
dc magnetic fields up to 35T. (b) \pxy\ vs pxx at various tem- 
peratures. The solid lines are fitting results using the scaling 
behavior \pj:y{H)\ = Apxx{H)^ . Inset of (b) shows the tem- 
perature dependence of PiT). 



samples, where /3 was found to decrease from 1.5 com- 
pared to 2 after irradiation,^^ and in HgBa2CaCu206+2; 
thin films with columnar defects, where /? changes from 
1.0 to 1.2 with increasing the fieldi^ Therefore, small 
values of /3 in S-09 (inset of Fig. 4(b)) may be connected 
with the strong pinning strength due to the considerably 
large concentration of defects in Fei+j,(Tei_a;Sa;)2i^ 



IV. CONCLUSION 

In summary, we investigated the resistive TAFF and 
flux flow Hall effect of Fci 14(1) (Tcq. 91(2)80. 09(2))z sin- 
gle crystal in high and stable magnetic fields up to 
35 T. TAFF behavior could be understood within the 
framework of modified Arrhenius relation assuming that 
the prefactor po/ is temperature-dependent while pc is 
temperature-independent, and U{T,H) = Uo{H){l — t)'^, 
q can be set as a free parameter for fitting not lim- 
ited to 1. There is a crossover from single- vortex pin- 
ning region to collective creep region for both field di- 
rection. Furthermore, q changes from 1 to 2 when mag- 
netic field is rotated along ab plane to c axis but it is 
not sensitive to the magnitude of magnetic field below 
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35T. Hall and longitudinal resistivity in mixed state in- 
dicate that there is no Hall sign change, as opposed to 
cuprate superconductors. We observed scaling behavior 
\Pxy{H) \ — Apxx{H)^ with the scaling exponent (3 about 
1, which may be due to strong pinning strength in con- 
siderably disordered in Fei+j,(Tei_a;Sa;)z system. 
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